Measuring spin correlations in optical lattices using superlattice potentials 
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We suggest two experimental methods for probing both short- and long-range spin correlations 
of atoms in optical lattices using superlattice potentials. The first method involves an adiabatic 
doubling of the periodicity of the underlying lattice to probe neighboring singlet (triplet) correlations 
for fermions (bosons) by the occupation of the new vibrational ground state. The second method 
utilizes a time-dependent superlattice potential to generate spin-dependent transport by any number 
of prescribed lattice sites, and probes correlations by the resulting number of doubly occupied sites. 
For experimentally relevant parameters, we demonstrate how both methods yield large signatures 
of antiferromagnetic (AF) correlations of strongly repulsive fermionic atoms in a single shot of the 
experiment. Lastly, we show how this method may also be applied to probe d-wave pairing, a 
possible ground state candidate for the doped repulsive Hubbard model. 

PACS numbers: 37.10.Jk, 73.21. Cd, 74.25.-q, 75.50.Ee 



The study of ultracold atoms in optical lattices has pro- 
duced several groundbreaking results including the obser- 
vation of quantum phase transitions to a Mott state [I], 
fermionic pairing [5] , and fermionization of bosons in one 
dimension [3]. A major, but presently unrealized goal, 
is to study quantum magnetism using atoms in optical 
lattices. For this purpose, it is important to have effi- 
cient methods to measure the atomic correlations. Possi- 
ble experimental probes include Bragg scattering [3] and 
quantum noise spectroscopy, where higher order corre- 
lation functions are measured by analyzing the fluctua- 
tions in time-of-flight [SHZj or light polarization exper- 
iments [8j. Time-of-flight experiments have been im- 
plemented experimentally to observe bosonic/fermionic 
bunching/anti-bunching and pairing 9J. The exper- 
imental signature for the spin correlations with this 
method is, however, very small and requires averaging 
over many shots of the experiment. Recent experiments 
have achieved single site resolution in probing optical lat- 
tice systems in two dimensions (2D) [TU1 E] , which could 
be used to probe correlations. 

In this Letter, we investigate two methods for measur- 
ing correlations in optical lattices of any dimensions us- 
ing superlattice potentials. Most importantly, these tech- 
niques allow for read-out of the correlations in the mean 
value of the relevant observable, and hence measure the 
correlations in a single shot of the experiment. The first 
method was realized in Ref. [T^] , and is sketched in Fig. 
[lja). A superlattice is slowly turned on while the original 
lattice is turned off thereby merging two neighboring lat- 
tice sites. As shown below, a measurement of the popula- 
tion of the new vibrational ground state directly provides 
the nearest-neighbor singlet population. This method is 
ideal for measuring short-range correlations typically ap- 
pearing at high temperatures T . Related methods were 
presented in [13l [M] allthough with much lower signals. 
To detect phase transitions, we introduce a new method 



to probe long-range correlations based on superlattice- 
induced spin-selective particle transfer by a controllable 
number of lattice sites, using the procedure illustrated 
in Fig l(c-e). A subsequent measurement of the fraction 
of doubly occupied sites Q] US] reveals the long-range 
spin correlations. As particular examples, we demon- 
strate how the method can detect AF and ci-wave pair 
correlations for two-component fermions in a lattice. 

We consider atoms with two internal levels |f) and 
\l) trapped in a cubic lattice of period a. It is suf- 
ficient to consider a single dimension x, assuming the 
motion in the other two dimensions to be frozen by 
large lattice potentials. Along the x-direction, the atoms 




FIG. 1: (Color online) (a) Adiabatically doubling the lattice 
wavelength forces the singlet (triplet) to enter the ground 
state for fermions (bosons). (b) Singlet probability for 
nearest-neighbors vs. T for the isotropic Heisenberg model, 
(c-e) A moving superlattice induces Landau-Zener transitions 
between neighboring lattice sites for spin-up atoms which have 
been transferred to the first excited states where the tunneling 
is higher than in the ground state. Atoms on even (odd) sites 
move right (left) following the superlattice minima (maxima). 
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are trapped by a time-dependent potential given by 
A{x,t) = A(t) sin 2 (7ra/a) [16] . In addition, we super- 
impose a superlattice potential of the form B(x,t) = 
B(t) sm 2 (Trx/2a — ttS/2). For measuring short-range cor- 
relations, two neighboring lattice sites are merged by 
adiabatically decreasing A(t) and ramping up B{t) with 
S = 1/2 (see Fig. [jja)). If this process is performed 
much faster than the tunneling time, but slower than the 
(much faster) oscillation period of the potential wells, one 
will adiabatically merge the two wells located at x = 2n 
and x = 2n + 1, n = 0, ±1, . . .. Two bosonic (fermionic) 
atoms in a singlet (triplet) state cannot occupy the vibra- 
tional ground state of the potential well due to the Pauli 
exclusion principle. Hence the singlet (triplet) state adia- 
batically connects to the first excited vibrational state of 
the two atom motional wavefunction. Thus, a measure- 
ment of the population in the vibrational ground state by 
expansion imaging reveals the probability for the atoms 
to be in the triplet (singlet) state. 

As an application of this technique, we consider a two- 
component Fermi gas in an optical lattice. For strong 
repulsion and half-filling, the gas is in the Mott phase at 
low T and is well described by the Heisenberg model 

ff=J^s r s ro . (1) 

(Lrn) 

Here, S; is the spin-1/2 operator for atoms at site I 
and (l,m) denotes neighboring pairs coupled by an ex- 
change interaction J. In three dimensions (3D), Eq.Q 
exhibits a phase transition to an AF ordered state at 
T = T/v — 0.945J, the observation of which is currently 
a main goal in the study of optical lattices. However, cur- 
rent experiments reveal only limited information on the 
spin state of the atoms and therefore on the temperature. 

To investigate signatures of AF correlations we con- 
sider initially the high-T regime, T» J, and perform a 
high T expansion of the density matrix 

where N is the number of atoms, and P* m (P[ m ) projects 
onto the singlet (triplet) state of atoms I and to. This 
expression shows that interactions increase the probabil- 
ity of neighboring pairs to be in a singlet state by 3 J/ AT. 
The probing method illustrated in Fig. [lja) is therefore 
ideal to measure the onset of AF correlations at high T. 

To investigate the behavior at lower T, we perform 
full Quantum Monte Carlo (QMC) simulations using the 
stochastic series expansion method [T7] with directed- 
loop updates [18] to calculate the correlation functions 
(si ■ Sj). This method is very efficient for Heisenberg 
models and gives accurate results for a wide range of 
T for large systems. Using the spin-correlations from 
the QMC simulations, we show in Fig. [jjb) the nearest- 



neighbor singlet probability vs. T. Evidently adiabati- 
cally merging wells is very efficient at high T. It does 
not, however, exhibit a clear signature for the phase 
transition at T/v, since it does not probe the onset of 
long-range correlations given by a non-vanishing value of 
= lim^oo C z m with C* = (sfsf +m ). 

This leads us to the second method which is to detect 
by selectively moving atoms m sites either to the left 
or right if they are in state |t), whereas ||) atoms re- 
main stationary. Assuming translational invariance, the 
probability for two atoms to be at the same site after 
this displacement is given by Pa = 1/4 — C^. A possible 
method to perform such a selective displacement based 
on spin selective optical lattices was proposed in Ref. [TS] 
and realized in Ref. [2U]. This procedure, however, re- 
quires lasers to be detuned less than the fine structure 
splitting of the atoms. For alkali fermionic atoms used 
in experiments this splitting is small and may cause sig- 
nificant heating [3T]. Our method does not suffer from 
this limitation, since it is based on optical superlattices 
which do not need to be spin dependent and therefore 
can have a large detuning. 

The central idea is to exploit that tunneling to neigh- 
boring sites is larger for the first excited state than for the 
vibrational ground state of each potential well. To dis- 
place only the |f) component, these atoms are excited to 
the first excited state, using e.g. a vibrational sideband 
of a Raman transition to an auxiliary state |o) (22] , while 
the 14.) state remains in the ground state. To selectively 
move |t) we then ramp up the superlattice B(x,t) with 
(5 = corresponding to potential minima at the even 
sites x = 2na of the underlying lattice. If we then let 
the superlattice potential move with S = vt/a, |f) atoms 
will be dragged along undergoing Landau-Zener transi- 
tions as illustrated in Fig. [TJc-e) : the vibrational eigen- 
energies at even and odd sites will anti-cross at 6 = 1/2, 
where the system basically consists of several two-level 
systems. Hence, the probability to cross diabatically P 
(remain on the same site) is given by the Landau-Zener 
formula, P = exp(— 2irg 2 /h\a\), where g is half the differ- 
ence between the energies at the avoided crossing, and a 
is the time-derivative of this energy difference [22] ■ Since 
the excited states in the neighboring wells are stronger 
coupled than the ground states (have larger g) the super- 
lattice velocity v s can be chosen so that the vibrational 
ground state remains stationary (P « 1), while the first 
excited state is moved (P « 0). Atoms in |f) at even 
(odd) sites thus move to the right (left) and the super- 
lattice acts as a spin-dependent "conveyor belt" . 

To demonstrate the feasibility of this proposal we 
model a system of 32 sites with lattice potential depths 
A = 40Er and B = 2QEr. Minimizing the sums of 
the errors (1 — P for the ground state and P for the ex- 
cited state) yields v s — 8.95 • 10~ 5 aE R /h. Since B/A 
is rather large, there are also level crossings between 
the ground state and the first excited state at different 



3 



sites. To ensure that these are diabatic (non-mixing) we 
move the superlattice with a high velocity v = aE^/h 
at these crossings. The velocity is decreased to the op- 
timal value only in an interval St > 0.02a/t> around the 
desired site-mixing anti-crossings. The total runtime is 
approximately T ~ m ■ 250h/Eji = m ■ 17 ms (m • 8 
ms) for the experimentally relevant case of 6 Li [2] ( 40 K 
[9]), with m being the number of moves we wish to per- 
form. In order to mimic an experimentally realistic situ- 
ation, we start by ramping up the superlattice potential 
B, and end the procedure by a similar ramp-down. The 
time-dependent Schrodinger equation is solved for each 
Bloch state of the superlattice unit cell by numerical in- 
tegration expanding in the 12 lowest eigenstates of the 
unperturbed lattice with B = 0. Transforming to a suit- 
able Wannier basis allows us to extract the transition 
matrices U(t) = e lHt = T,ij„ m U ijnm(t)at n a im) where 

a\ n is the creation operator for the nth vibrational mode 
Wannier function located at site i. The time evolution 
operator U(t) is quadratic in a and since we assume 
that the interaction effects can be tuned to zero using 
a Feshbach resonance. We utilize Floquets' theorem for 
periodic Hamiltonians to extract the time-evolution. 

Fig. [2|a) shows the norm square \Wi\ 2 of the calcu- 
lated Wannier functions after m = 1, . . . , 8 moves, for an 
atom starting in the first vibrational state at an even site. 
For this state the procedure works very well. The inset 
in Fig. |2ja) confirms that the ground state is immune 
to the moving superlattice. Figure [2jb) shows the corre- 
sponding odd site Wannier function which moves to the 
left. The inset in Fig. |2jb) gives the error, i.e. the prob- 
ability to be moved anything but m sites from the initial 
position or leak to other vibrational states. Note that 
for the chosen parameters, the error is higher for the left 
moving atoms as compared to the right moving atoms, 
even though the two-site Landau-Zener problem is sym- 
metric with respect to the left and right moving atoms. 
This is because the left moving atoms follow the maxima 
of the superlattice potential B(x,t) making them more 
susceptible for tunneling to states two sites away as seen 
by the growing "sidepeaks" at m± 2 in Fig. [2] This res- 
onant tunneling can be quenched by using e.g. a larger 
B at the cost, however, of a smaller signal because of 
enhanced transitions to higher vibrational modes. 

Assuming the (de)excitation process for |f) to the first 
vibrational state is perfect, the probability Pd(i) for site 
i to have two atoms in the ground state n = is 

PS) = E fel 2 l^ooi 3 (i + 1 (s'i - st) - (spi)) 

~ X, u *jii u ii^ u mo u iW {SjSj) . (3) 

Note that a measurement probing all sites leads to a P#- 
signal proportional to the system size N. Due to the dif- 
ference between states originating in odd and even sites, 
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FIG. 2: (Color online) (a) Probability distribution | Wi\ 2 of the 
first vibrational excitation at even sites for moves m = 1, . . . 8. 
The curves are offset for clarity. The inset confirms that the 
vibrational ground states remains stationary, (b) Probability 
distribution of the first excited states on odd sites. The inset 
shows the error probability for the three states in (a) and (b) . 



we consider the averaged (Pd(i)) = [Pd(i) + Pd{i + l)]/2- 
In Fig. [3^a-c), we show (Pd(i)) function of T for 

different moves m for the Heisenberg model in ID, 2D 
and 3D. The curves are obtained from Eq.(|3| by combin- 
ing the calculated transition matrices V^m for a given 
m with (§; • Sj) obtained from QMC simulations for the 
isotropic Heisenberg model. The increase in spin corre- 
lations with decreasing T clearly shows up as an increase 
in (Pd(i)) for m odd and a decrease for m even. The 
correlation decreases with increasing site separation m, 
and Fig. [3^c)-(d) show how the onset of long-range order 
for T < T/v for a 3D system results in a dramatic change 
in (Pd(i)} for large m as a direct consequence of the spin 
ordering in the AF state. Note that (Pd(i)) does not con- 
verge to the uncorrelated value 1/4 when T/J — > oo since 
the transfer procedure leads to atom loss into higher vi- 
brational bands, which we do not include. As mentioned 
above, the pile-up of weight at the sidepeaks m±2 in Fig. 
[2] can be minimized using other parameters, allowing a 
more direct comparison between (Pd(i)) and C^. Note, 
however, that this contribution to the "error" does not 
contribute negatively to (Pd(i)) in the AF state because 
m ± 2 belong to the same sublattice as m. 

Fig. ^d) shows (Pd(i)) vs. T in 3D in the presence 
of a small exchange anisotropy. This case reveals the 
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FIG. 3: (Color online) Average probability (Pd(i)) for a site 
to be doubly occupied vs. T after moves m = 1, . . . , 8 for the 
isotropic ID (a), 2D (b), and 3D (c) Heisenberg model, (d) 
same as (c) with the addition of a small exchange anisotropy 
H z =0.01,7 J2 ([m} sts z m . 



strongest signature of the AF phase transition with a 
pronounced difference for large number of moves, since 
we force the symmetry breaking to be along the direction 
that we probe. By contrast, for the ID and 2D results in 
Fig. [3];a)-(b) the curves only "peel off" the (P d (i)) = 1/4 
line for decreasing T, since there is no phase transition. 
Using this method, it is thus possible to probe long-range 
correlations, and obtain large experimental signals in a 
single shot of the experiment. 

Above we focused on the detection of AF correlations, 
but the method can be used to probe any correlations 
in optical lattices. A particularly interesting case is the 
possibility of d-wave pair correlations within the phase 
diagram of the repulsive Hubbard model relevant to the 
field of high-temperature superconductors. Such correla- 
tions are encoded in the d x 2_ y 2-w&ve correlation function 

Cfj = (A- Aj)> where d\ = a\i, & i+^~ & iA+m with * (w) 
unit vectors along the x (y) direction. To measure C D , 
we assume that an initial operation is performed which 
swaps |t) atoms on neighboring sites along the x-axis. 
Next a pulse of duration St is applied, e.g. a photoassoci- 
ation pulse, merging two atoms into a molecule if they are 
on the same site. This process is described by the Hamil- 
tonian H = ^ J2i b\a^cii^ + FT. a, where f2 denotes the 
strength of the interaction and bi is the bosonic annihi- 
lation operator for the molecular field. After the applied 
pulse, unassociated atoms are transferred back along the 
x-direction while molecules remain at rest. Neighbor- 
ing |t) atoms are then swapped along the y-direction, 
and another association pulse is applied with the oppo- 
site phase £1 — > —CI. Finally, unassociated atoms are 



moved back along y. Solving for the molecular field af- 
ter the pulse we find b° ut = bf - iVt8tD l /2. Thus, the 
probability to find a molecule at a particular site is di- 
rectly related to the probability of finding a d-wave pair 
on this and the neighboring sites. The total number of 
molecules N b is obtained to be N b = Sl 2 (5t 2 C t ?/ 4 = 
fl 2 St 2 X; 4 (2n lt n 4 +4|A i | 2 )/4, where n l(J = (a] CT a ier ), and 
A, = {cLi+Sxhl) = -{o-i+ino-il) denote the density and 
d-wave pairing gap, respectively. An expansion image 
of the molecular cloud should therefore reveal <i-wave 
pairing as a peak at zero momentum [5J [7J. The pres- 
ence/absence of this peak at T > T c can also help eluci- 
date whether the pseudo-gap phase of high-T c cuprates is 
caused by superconducting fluctuations or rather a hid- 
den (non-superconducting) order. 

In summary, we have demonstrated how to use super- 
lattice potentials to probe both short-range and long- 
range correlations in optical lattices. Specifically, we 
showed how AF and d-wave superfluid correlations can be 
measured in a single shot of the experiment by this pro- 
cedure. The method is, however, applicable as a probe of 
any correlations, and should be highly useful for future 
studies of quantum many-body systems. 
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